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ABSTRACT 


The  research  conducted  under  this  contract  concerned 
eight  general  areas  of  investigation: 

1.  Establishment  of  a  network  of  seismic  stations 
in  the  central  U.S.  suitable  for  the  study  of  long 
period  surface  waves  and  other  specialized  studies, 

2.  Phase  and  group  velocity  studies  of  surface 
waves.  These  include  phase  velocity*“de terminations 
of  both  Love  and  Rayleigh  waves,  inversion  of  the 
velocity  data,  group  velocities  of  very  short  period 
Rayleigh  waves,  and  group  velocities  along  paths 
which  cross  the  Caribbean  Sea, 

3.,  Reduction  of  ground  motion  to  motion  in  the  in¬ 
cident  wave.  Theoretical  work  by  Haskell  has  been 
verified  from  observational  data  and  further  developed 
to  study  angles  of  Incidence  of  P  waves  and  polariza¬ 
tion  of  S  waves  which  traverse  a  multilayered  crust. 
Crustal  transfer  functions  for  various  crustal  models 
have  been  determined  and  used  to  study  the  effect  on 
ground  motion  of  changes  in  crustal  parameters. 


4.  Focal  mechanism  determinations.  Radiation 
patterns  of  P  and  S  waves  are  found  to  be  characterie- 
tic  of  given  earthquake  sequences  and,  to  some  degree, 
of  given  earthquake  regions. 


5.  Focal  depth.  The  Pg-Pn  vs  sg~pg  raethcd  has 
been  found  of  qualitative  help  in  studying  focal 
depth  of  shocks  in  SE  Missouri,  The  transfer  func¬ 
tion  for  a  source  in  a  layered  medium  has  further 
shown  that  at  shallow  depths  tne  source-crust  system 
acts  as  a  rejection  filter  for  the  low  frequency 
portion  of  the  spectra  of  body  waves. 


6.  Inversion  of  body  wave  data.  The  transfer 
function  of  TW/TU  (ratio  of  vertical  to  horizontal 
ground  motion)  has  been  calculated  as  a  function  of 
a  dimensionless  frequency  for  representative  one- 
and  two-3 ayer  crystal  models.  Spectra  of  observed 
P  waves  are  then  compared  to  the  theoretical  curves 
to  determine  crustal  structure. 


7.  Effect  of  low  velocity  channels  on  travel 
tii-:es.  Parameters  of  a  low  velocity  channel  have 
been  systematically  varied  to  examine  the  circum¬ 
stance  under  which  the  channel  produces  or  does  not 
produce  a  shadow  zone,  BILBY  and  GNOME  data  were 


then  cx^lned  to  3tud.y  lateral  variations  in  upper 
mantle  velocities  throughout  the  central,  eastern 
and  western  U.  S. 

±  Lonal  crustal  structure,  seismicity,  Pn  veloc 
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Introduction 


The  research  conducted  under  this  contract  concerned 
eight  general  areas  of  investigation: 

1.  The  establishment  and  operation  of  a  quadri¬ 
lateral  network  of  seismograph  stations  In  the 
central  United  States. 

2.  Phase  and  group  velocity  studies  of  seismic 
surface  waves . 

3.  Reduction  of  ground  motion  to  motion  in  body 
waves  incident  at  the  base  of  the  crust. 

4.  Focal  mechanism  studies. 

5.  Focal  depth  of  seismic  sources. 

6.  Inversion  of  body  wave  spectra  to  determine 
crustal  structure. 

7.  Effect  of  low  velocity  channels  on  travel 
times  of  seismic  waves. 


B.  Local  crustal  structure,  seismicity,  Pn  velocity. 
The  results  which  have  been  forthcoming  from  these  investi¬ 
gations  have,  for  the  most  part,  been  published  in  a  series 
of  papers.  A  list  of  these  papers  may  be  found  at  the  end 
of  this  report.  In  the  body  of  the  report  we  shall  first 
give  a  brief  discussion  of  the  results  of  the  research 
under  each  of  the  areas  of  investigation  mentioned  above. 
This  discussion  will  be  followed  by  the  abstracts  of  the 


papers 


1.  Discussion  of  Results. 


1.1  Network  of  Seismograph  Stations.  A  quadrilateral  net¬ 
work  of  seismic  stations ,  similar  in  instrumentation  to  the 
V/WSS  system  was  established  at  Manhattan,  Kansas,  Rolla, 
Missouri,  Dubuque,  Iowa,  and  Bloomington,  Indiana.  A  seizes 
of  calibration,  experiments  was  performed  in  order  to  select 
the  Instrumental  constants  most  appropriate  for  the  record¬ 
ing  of  S  waves  and  mantle  Rayleigh  waves,  that  is,  waves  in 
the  period  range  of  10  seconds  to  100  seconds  and  longer. 

The  network  has  been  essential  for  the  supply  of  data  for 
the  long-period  P,  S,  and  surface  wave  studies  described 
below. 

In  addition  to  contract-related  work,  copies  of  records 
and  other  data  have  routinely  been  supplied  to  the  USCGS  and 
to  other  VELA-Unlform  users  for  NTS  and  events  of  similar 
type.  A  phenomenon  which  has  been  noted,  but  thus  far  not 
Investigated  further,  is  an  early  arrival,  prominent  on  the 
short^period  vertical  record  at  Manhattan,  for  NTS  events. 
This  arrival  Is  about  4  seconds  earlier  than  the  expected 
time  of  arrival  of  Pn;  the  latter  phase  is  also  a  prominent 
arrival. 

Two  long-playing  short-period  seismic  stations  nave 
been  installed  in  trailers  for  temporary  use  in  the  SE 
Missouri  seismic  zone,  and  for  the  recording  of  other 
events  of  interest. 
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1*2  Phase  -tnd  Group  Velocity  Studies.  Our  primary  studies 
concerned  phase  velocities  in  a  mid-continental  region. 

These  studies  have  been  supplemented  by  certain  other  sur¬ 
face  wave  investigations. 

In  the  central  United  States,  over  the  area  of  the 
St.  Louis  University  network,  phase  velocities  of  both  Ray- 
Lgh  and  Love  waves  have  been  measured  for  the  period  range 
5  to  80  seconds.  These  data  are  of  significance  for  two 
reasons.  First,  while  the  Love  wave  phase  velocities  in 
the  period  range  considered  were  found  to  be  comparable  to 
those  determined  by  Brun~  and  Dorman  for  the  Canadian  Shield 
the  Rayleigh  wave  phase  velocities  are  notably  lower  (3$) 
in  the  central  U.S.  than  those  in  the  Canadian  Shield.  Sec¬ 
ond,  an  inversion  program  designed  to  determine  a  crust- 
upper  mantle  model  which  would  fit  both  the  Love  and  the 
Rayleigh  wave  phase  velocity  data  was  unsuccessful  in  achiev¬ 
ing  a  simultaneous  fit  to  the  two  kinds  of  data.  Aniso¬ 
tropy  seems  to  be  required  in  the  upper  mantle.  This  is  in 
contrast  to  the  situation  which  prevails  in  the  Canadian 
Shield.  These  results  indicate  a  need  for  further  work  in 
order  to  investigate  r.)  the  transition  between  the  Canadian 
Shield  structure  and  the  central  U.S.,  b)  the  dimensions 
of  areas  ever  which  phase  velocities  are  constant,  c)  the 
evidence  and  degree  of  anisotropy  in  the  upper  mantle. 

A  study  of  short-period  (q.2  to  2.0  sec.)  Rayleigh 
wave  dispersion  in  signals  from  quarry  and  strip  mine  blasts 
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and  small  jocal  earthquakes  has  shown  the  sensitivity  of 
Rayleigh  waves  of  this  period  range  to  changes  in  the  sedi¬ 
mentary  station.  This  suggests  their  use  as  a  tool  in 
shallow  crustal  x*econnaissance . 

Group  and  phase  velocities  along  57  paths  from  central 
America  co  central  U.S.  were  also  used,  applying  establisnea 
methods  of  interpretation,  in  determining  differences  in 
structure  between  the  western  and  eastern  Gulf  of  Mexico 
and  the  Caribbean  Sea . 

1.3  Reduction  of  Ground  Motion  to  Motion  In  Incident  Wave . 
Any  study  of  body  wave  amplitudes  which  is  concerned  with 
the  source  mechanism  or  the  velocity  structure  of  the  earth's 
interior  must  take  account  of  the  amplitude  and  phase  distor¬ 
tion  caused  by  crustal  layering  and  the  free  surface  at  the 
recording  site,  as  well  as  that  caused  by  the  seismograph. 

A  knowledge  of  the  response  characteristics  of  the  seismo¬ 
graph  is  sufficient  to  compensate  for  the  latte ?,  i.e.,  to 
reduce  the  trace  motion  on  the  seismogram  to  ground  motion 
at  the  recording  site.  The  other  aspect  of  the  problem,  to 
relate  ground  motion  tc  earth  motion  beneath  the  crust  near 
the  recording  site,  has  occupied  much  of  cur  attention. 

1.31  P  Waves .  The  relation  of  the  ground  motion  to  the 
P  wave  particle  motion  beneath  the  seismograph  station  has 
attracted  the  attention  of  many  renowned  seismologists,  such 
as  Galitsin,  Byerly  and  his  students,  Gutenlerg,  and  Hodgson. 
The  reasons  for  their  interest  in  this  relationship  were 


quite  different,  and  Include  suen  basic  problems  in  seis¬ 
mology  as  the  detection  and  location  of  second-order  veloc¬ 


ity  discontinuities  in  the  earth's  mantle,  the  calculatjon 
of  f-  jrtenried  i 

*  -  —  -**  *i«.v»iuiuoiu  dUUUlvOj  O I  ItT  &  is<lm 

‘  correction  for  magnitude  studies,  and  the  effect  of 
unconsolidated  materials  on  earthquake  intensity  and  pos¬ 
sible  structural  damage.  To  these  may  be  added  a  problem 
which  hap  received  much  attention  in  the  VELA-Uniform  pro¬ 
gram,  to  determine  station  sites  which  have  a  large  signel- 
to-noise  ratio  for  the  P  wave  and  to  determine  how  this 


ratio  changes  at  an  individual  location  with  epicentral  dis¬ 
tance  and  azimuth. 


One  result  of  the  research  undertaken  in  this  contract 
was  to  verify  experimentally  Haskell's  conclusions,  based  on 
a  theoretical  study,  that  the  P-wave  ground  motion  depends 
on  the  crustal  layering  at  the  seismograph,  station „  This 
was  accomplished  by  calculating  the  apparent  angle  of  inci¬ 
dence  of  the  P  wave  at  the  seismograph  stations  of  our  net¬ 
work  for  teleseismlc  events.  In  addition,  the  dominant 
period  of  the  first  half-cycle  of  the  P  vrave  ground  motion 
was  measured,  in  order  to  determine  the  frequency  dependence 
of  the  relationship  between  the  apparent  angle  of  incidence 
and  the  epicentral  distance.  The  observational  data  com¬ 
pared  well  with  values  calculated  from  Haskell's  paper, 
and  demonstrated  the  necessity  of  including  the  frequence- 
dependence  in  such  studies.  In  audition,  the  study  explained 


the  earlier  and  seemingly  paradoxical  observations  of  Mei 
and  Fyerly,  and  of  Nuttli  and  Whitmore,  that  in  order  to 
satisfy  the  short-period  data  with  a  half-space  model  of 
the  earth  it  was  necessary  to  assume  a  velocity  of  about 
8  km/sec.,  whereas  Papasachos  required  a  velocity  of 
6.4  km/sec.  to  satisfy  P-wave  data  with  periods  of  7  to  15 
se.onds,  Thus,  both  Haskell's  theoretical  calculations 


ment  which  appeared  so  often  in  the  seismological  literature 
to  th  point  of  almost  being  axiomatic,  that  the  effective 
half-space  velocity  for  the  short-period  waves  is  that  of 
the  upper-crustal  layers  whereas  that  for  the  long-period 
waves  is  the  upper-mantle  velocity. 

Hannon  extended  Haskell ' s  work  by  considering  more 
complex  crustal  models.  He  selected  models  which  demon¬ 


strated  the  effect  of  changes  in  crustal  thickness,  the 
presence  of  l  low- velocity  layer  of  variable  thickness,  such 


as  alluvium,  at  the  outer  surface,  and  the  presence  of  a 


low- velocity  layer  ir.  the  interior  of  the  crust.  Several 


Important  conclusions  follow  from  this  study.  These  are: 
l)  The  presence  of  a  low-velocity  material  at  the  surface 
causes  the  ground  motion  to  be  amplified.  Thus  one  can  ex 


pect  a  P-wave  signal  to  have  a  greater  amplitude  at  a  re¬ 
cording  site  set  up  on  lower-velocity  sedimentary  rocks 
than  at  one  on  higher-velocity  igneous  rocks.  2)  The  trans 


mission  coefficients  for  P  waves  with  periods  greater  than 
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abo,,J-  ,.w<  seconds  principally  depend  on  the  total  crustal 
thickness,  the  average  velocity  of  the  crustal  layers,  and 
the  contrast  between  the  upper  mantle  and  the  average 
crustal  velocity.  Thus, in  doing  the  inverse  problem  of 
determining  a  crustal  model  from  the  amplitude  spectrum  of 
the  observed  ground  motion,  it  is  only  possible  to  deter¬ 
mine  gross  crustal  properties,  such  as  thickn-ss  and  aver¬ 
age  velocity,  and  not  the  fine  structure,  such  as  the  dif¬ 
ferentiation  between  velocity  gradients  and  velocity  dis¬ 
continuities,  unless  one  uses  frequencies  of  one-half  cycle 
per  second  anu  higher.  3)  The  presence  of  an  internal  low- 
velocity  layer  ir  the  crust  has  an  effect  on  the  transmis¬ 
sion  coefficients  only  in  the  higher-frequency  range. 

Studies  of  the  ratio  of  the  amplitude  of  the  PP  and 
P  wave  ground  motion,  which  are  an  extension  of  earlier  work 
of  Gutenberg  and  Richter,  indicate  no  essential  difference 
in  the  crustal  structure  of  the  Atlantic  and  Pacific  Oceans. 
Both  the  amplitude  ratios  and  the  existence  of  two  separate 
PP  waves,  one  reflected  from  the  base  of  the  crust  and  the 
other  from  the  free  surface,  are  useful  in  identifying  a 
region  at  which  the  reflections  take  place  as  being  con¬ 
tinental  in  character,  as  opposed  to  oceanic. 

Long-period  P-wave  ground  motion  (periods  of  6  or  7 
seconds  and  greater)  from  tele seisms  has  been  found  to  be 
essentially  the  same,  both  in  amplitude  and  in  wave  form, 
for  all  stations  of  our  central  United  States  network  for  a 


given  earthquake.  This  indicates  that  the  crustal  struc¬ 
ture  is  almost  uniform  throughout  this  area,  and  also  that 
the  sedimentary  layers,  which  vary  considerably  in  thick¬ 
ness  and  in  rock  type,  have  no  discernible  effect  on  the 
long -period  P-wave  ground  motion. 

1.32  S  Waves.  Studies  of  the  amplitude  and  phase 
relationship  of  the  components  of  the  S-wave  ground  motion 
were  pursued  for  two  purposes.  The  first  concerns  the  S- 
wave  polarization  angle,  which  has  been  used  with  much  suc¬ 
cess  in  source  mechanism  studies.  The  second  relates  to 
the  velocity-depth  variation  of  S  in  the  upper  mantle,  in¬ 
sofar  as  it  affects  the  amplitude  and  travel  time  of  S. 

Our  observational  studies  of  the  S-wave  ground  motion 
shewed  that  this  ground  motion  was  linear  for  eplcentral 
distances  larger  than  about  forty  to  forty- five  degrees, 
and  elliptical  and  quite  complicated  for  lesser  distances. 

It  is  important  tc  fix  the  distance  at  which  the  ground 
motion  changes  from  non-linear  to  linear,  for  it  is  only 
at  the  greater  distances  that  the  polarization  angle  of  the 
S  wave  is  approximately  equal  to  the  angle  between  the  direc¬ 
tion  of  the  horizontal  component  of  the  S-wave  ground  motion 
and  the  great-circle  path  at  the  seismograph  station. 

Haskell's  theoretical  work,  which  considered  the  in¬ 
fluence  of  crustal  layering  at  the  seismograph  station, 
showed  that  the  above-mentioned  observations  agreed  with 
his  calculated  values  for  a  continental-crust  earth  model. 
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He  also  noted  that  ti a  amplitude  and/or  phase  of  the  SV- 
component  of  the  ground  motion  varied  rapidly  with  wave 
frequency  and  angle  of  incidence  at  the  lesser  eplcentral 
distances.  Although  this  would  indicate  that  the  polariza¬ 
tion  angle  cannot  be  obtained  from  seismograms  at  distances 
less  than  forty-five  degrees ,  hi3  analysis  suggested  a 
method  which  enabled  us  to  solve  the  problem.  By  a  slight 
extension  of  Haskell's  calculations,  we  have  been  able  to 
construct  a  set  of  charts  that  permit  one  to  obtain  the 
polarization  angle  from  the  ground  motion  for  distances  as 
small  as  twenty  degrees.  Stauder  and  Bollinger,  in  their 
study  of  the  source  mechanisms  of  the  large -magnitude  earth¬ 
quakes  of  1962  and  1963,  conclude  that  polarization  angles 
calculated  by  using  these  charts  for  distances  less  than 
forty-five  degrees  are  in  conformity  with  polarization 
angles  at  the  larger  distances  which  were  obtained  in  the 
conventional  manner. 

The  study  of  amplitude  variations  and  travel  times  of 
S- waves  which  have  traveled  a  large  part  of  their  path  in 
the  upper  mantle  is  difficult,  because  these  waves  arrive 
at  the  eplcentral  distances  for  which  the  3-wave  ground  mo¬ 
tion  is  extremely  complicated.  A  method  of  determining  the 
onset  time  of  S  has  been  developed,  and  is  presently  being 
applied  to  the  data  obtained  at  the  LHSM  stations  for  about 
eight  earthquakes. 


1C 
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*  Early  in  the  TELA-Uniform  program  a  suggestion  was 

made  that  amplitude  ratios  such  as  P/SH  or  SV/SH  could  pro- 

t 

vide  a  diagnostic  tool  for  distinguishing  between  explosions 
and  earthquakes,  because  it  was  expected  that  explosions 
would  produce  little  SH  motion.  Even  if  this  premise  were 
true,  one  cannot  apply  the  method  without  reducing  the 
ground  motion  to  particle  motion  within  the  earth.  For 
example,  the  crustal  layers  and  the  free  surface  of  the 

I 

earth  have  such  a  large  effect  upon  the  S-wave  motion  that 
a  wave  whose  particle  motion  is  almost  exclusively  SV  be- 
neath  the  crust  at  the  station  can  produce  a  ground  motion 
which  is  almost  completely  SH  in  character.  The  conclusion 
which  follows  is  that  ratios  of  P/SH  or  SV/SH  of  the  ground 
motion,  i.e.  motion  corrected  for  the  instrument  response 
characteristics  but  not  for  the  crustal  layering  effects, 
are  of  little  value  and  may  even  be  misleading  when  used 
as  diagnostic  techniques  to  classify  an  event  as  an  explo¬ 
sion  or  an  earthquake. 

1.4  Focal  Mechanism  Studies .  Our  focal  mechanism  studies 
have  progressed  by  the  application  of  the  theory  of  point 
sources  to  the  interpretation  of  observational  data  of  the 
P  and  S  waves  from  teleseisms.  Work  originally  begun  under 
Contract  AF  19(604)-8054  has  been  incorporated  into  this 
effort.  We  have  also  examined  groups  of  fault  plane  solu- 
.  tions  by  many  authors  for  earthquakes  around  the  borders 

of  the  Pacific  Ocean. 
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We  have  demonstrated  that  S  wave  polarization  infcrma- 
*  tlon  can  be  used  as  a  powerful  tool  in  the  determination  of 
focal  iuechanisms.  While  S  wave  data,  in  our  opinion,  are 
most  advantageously  used  as  an  adjunct  to  the  ?  wave  first 
motion  technique,  we  have  also  developed  separate  method3, 
using  S  wave  data  alone,  A  computer  program  has  been  de¬ 
veloped  to  this  end,  and  has  been  used  effectively  in  study¬ 
ing  a  number  of  earthquakes. 

Certain  radiation  patterns  of  P  and  S  waves  have  been 

#  t 

found  to  be  cliarabterlstic  of  a  given  earthquake  sequence 
or  of  earthquakes  in  a  given  region.  Of  particular  interest 
to  detection  problems,  the  occurrence  of  compressional  first 
motion  of  P  waves  for  earthquakes  as  recorded  in  the  third 
zone,  together  with  SV  polarization  of  the  S  waves,  is  not 
a  rare • occurrence .  Earthquakes  of  the  Kamchatka  area  in  pap 
ticular  are  commonly  of  this  nature. 

The  limitations  inherent  in  methods  of  focal  mechanism 
determination  from  the  data  of  body  waves  have  been  examined 
Chief  among  these  are  the  limitation  in  geographic  distribu¬ 
tion  of  stations  and  the  small  dimensions  the  seismic 
window  on  the  focal  sphere  through  which  it  is  possible  to 
see  the  earthquake  focus  from  data  of  body  waves  recorded 
at  Pn  distances  and  greater. 

Examination  of  multiple  solutions,  i.e.,  by  several 
different  authors  for  one  and  the  same  earthquake,  as  also 
a  comparison  of  groups  of  fau?.t  plane  solutions  for  a  given 


region  evidence  poor  agreement  and  systematic  differences  in 
solutions  by  different  authors.  This  evidence  requires  that 
caution  should  be  exercised  in  drawing  conclusions  from  the 
existing  accumulation  of  published  fault  plane  solutions. 

Progress  has  been  made  in  the  use  of  S  wave  polariza¬ 
tion  data  for  smaller  earthquakes  at  shorter  epicentral  dis¬ 
tances.  Use  of  phase  equalization  of  surface  waves  was 
attempted  for  a  magnitude  5  earthquake  of  southeast  Missouri 
but  the  results  were  not  conclusive. 

1.5  Focal  Depth.  Two  approaches  have  been  tried  to  the 
focal  depth  problem. 

By  taking  advantage  of  a  known  surface  event,  the  BILLI- 
KEN  shot  of  June  28,  1963,  the  Pg-Pn  vs  Sg-Pg  method  was  ap¬ 
plied  to  two  magnitude  5  earthquakes  of  southeast  Missouri, 

An  earthquake  of  2  February  1962  was  found  to  have  a  shallow 
focal  depth;  one  of  3  March  1963  to  have  a  focus  deep  12c  the 
crust.  It  was  also  noted  that  the  seismogram  of  the  latter 
carried  considerably  more  power  in  higher  frequency  portion 
of  the  record  than  did  the  former. 

Secondly,  applying  to  body  wave  generation  a  souree- 
c~  ust  model  set  up  by  Harkrider  for  surface  waves,  we  have 
obtained  the  transfer  function  for  a  source  in  a  layered 
system  overlying  a  half-space  as  a  function  of  source  depth. 

It  Is  found  that  for  an  explosive  point  source  located  at 
shallow  depths  the  source-crust  system  acts  as  a  rejection 
filter  for  the  low  frequency  portion  of  the  spectra  of  body 
waves.  This  conclusion  seems  to  point  to  a  most  promising 
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technique  of  focal  depth  determination  and  indicates  the 
advisability  of  further  work. 

1.6  Inversion  of  Body  Wave  Spectra.  Two  studies  have  con¬ 
tributed  to  the  problem  of  inferring  crustal  structures 
under  a  recording  site  from  the  spectx^a  of  P  waves  recorded 
at  the  site. 

In  a  preparatory  study  Hannon  (Technical  Report  3, 

AFCRL  64-614)  investigated  the  transfer  function  of  a  lay¬ 
ered  crust  and  examined  the  effect  of  varying  the  parameters 
of  the  crust  model.  He  obtained  synthetic  seismograms  for 
P  waves  recorded  at  the  surface  of  various  crustal  models. 

He  remarked  the  differences  of  the  seismograms,  but  did  not 
undertake  the  Inverse  problem. 

Following  the  lead  of  Hannon,  Fernandez  Investigated 
the  possibility  of  a  Least  squares  Inversion  whereby  a 
crustal  model  might  be  fitted  to  observed  spectra  of  P  waves. 
He  found  the  partial  derivatives  with  respect  to  velocity, 
layer  thickness,  etc.,  of  the  transfer  function.-  when  ex¬ 
pressed  as  a  function  of  frequency,  to  be  so  large  and  to 
vary  so  erratically  as  to  make  an  inversion  program  not 
feasible.  By  expressing  the  transfer  function  In  terms  of  a 
dimensionless  parameter,  however,  a  much  more  regular  be¬ 
havior  was  obtained.  While  he  did  not  develoo  a  least 
squares  Inversion  program,  Fernandez  has  calculated  a  set 
of  master  curves  for  the  transfer  functions  of  one  and  two- 
layer  crusts.  By  overlaying  observed  spectra  (or  the  ratio 
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TW/TU  of  the  spectrum  of  the  vertical  component  to  th*t  of 
the  horizontal)  on  the  master  curves,  a  best  fitting  crust 
is  found  by  a  curve -matching  technique.  Using  long-period 
records  of  the  Saint  Louis  University  net,  the  technique  is 
remarkably  successful  in  determining  the  average  features 
of  the  crust.  The  method  indicates  that  in  the  central  U.S. 
the  average  crustal  velocity  is  6.6  km/sec  and  the  crustal 
thickness  is  42  km.  Finer  details  of  crustal  structure  may 
be  determinable  using  the  shorter  periodpportion  of  the  body 
wave  spectra.  This  material  is  being  submitted  to  AFCRL  as 
Scientific  Report  Number  13. 

1 . 7  Effect  of  Low- Velocity  Channel  on  Travel  Times  of 
Body  Tvaves. 

From  the  pioneering  work  of  Gutenberg  and  others  it 
is  well-known  that  low-velocity  channels  in  the  earth  can 
produce  shadow  zones.  The  best  example  of  this  is  the  shadow 
zone  for  the  P  wave  resulting  from  the  earth’s  core.  There 
is  both  observational  evidence  and  geometrical  proof  obtained 
by  constructing  ray  paths  that  the  existence  of  a  region  in 
the  earth  for  which  r/v  increases  with  depth  can  produce  a 
shadow  zone  if  the  hypocenter  is  located  in  or  above  the  low- 
velocity  -channel.  But  too  little  attention  had  been  given  to 
the  converse  problem,  i.e.,  does  the  presence  of  a  low- 
velocity  channel  always  produce  a  shadow  zone.  A  survey  of 
the  literature  suggested  that,  except  for  the  work  of  Leh¬ 
mann  and  Nishlraura,  the  converse  was  generally  accepted  as 
being  true. 


The  problem  does  not  alio*  a  simple  analytical  solu¬ 
tion,  because  the  expressions  for  the  travel  time  and  epl- 
central  distance  of  a  ray  traveling  through  a  earth  model 
which  traverses  a  region  in  which  r/v  Increases  with  dem 

*  jr  — 

Involve  non-convergent  Integrals .  But  from  purely  physical 
reasoning  we  can  argue  that  there  must  be  a  solution,  be¬ 
cause  the  rays  encountering  such  a  low-velocity  channel  do 
arrive  at  the  earth «s  surface  in  a  finite  amount  of  time. 

The  method  we  used  to  solve  this  problem  employed  ray  track¬ 
ing  which  was  programmed  for  a  digital  computer. 

A  systematic  study  was  made  of  the  various  parameters 
of  the  low-velocity  channel  which  can  affect  the  time-distance 
curve.  The  most  significant  result  was  tha.  a  sufficiently 
large  velocity  gradient  (increase  of  velocity  with  depth) 
will  not  only  produce  no  shadow  zone,  but  will  give  an  over¬ 
lapping  of  the  segments  of  the  curve  so  that  it  is  possible 
to  have  three  wave  arrivals  at  a  given  distance,  a  result 
which  is  qualitatively  similar  to  that  produced  by  an  earth 
model  with  no  low- velocity  channel  but  rather  an  abrupt  in¬ 
crease  in  the  velocity  or  velocity  gradient  with  depth. 

Attention  was  also  given  to  the  inverse  oroblem,  namely 
to  determine  the  velocity-depth  relation  from  P-wave  travel 
times.  The  data  used  were  t he  P  times  in  North  America  from 
the  nuclear  events  BILBY  and  mm.  Although  there  is  some 
ambiguity  in  the  details  of  the  velocity-depth  relation 
obtained  in  this  manner,  it  can  definitely  be  stated  that 
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there  are  lateral  variations  in  the  upper-itantle  P-wave  veloc¬ 
ities.  The  central  and  eastern  parts  of  »;he  United  States 
have  in  general  higher  Pn  velocities  than  the  western  part. 
There  is  no  strong  evidence  for  the  existence  c.f  a  lew** 
velocity  channel  for  P  in  the  central  and  eastern  United 
States.  If  one  exists  it  is  relatively  deep*  not  very  thick, 
and  entails  only  a  smalx  velocity  decrease.  The  upper-mantle 
P-wave  velocity  structure  in  tnt  western  United  States^  how¬ 
ever*  is  more  complex  and  suggests  large  lateral  variations. 

The  low-velocity  channel  for  P  is  shallower  and  more  pro¬ 
nounced*  and  can  produce  shadow  zones  for  P  such  as  those 
postulated  by  Gutenberg  for  southern  California. 

1,8  Local  Crustal  Structure*  Seismicity*  Pn  Velocity.  Several 
different  studies  fall  under  this  heading. 

*'T 

A  continuing  study  of  he  seismicity  of  the  southeast 
Missouri  area  has  leu  to  new  understanding  of  this  main  seib- 
mic  region.  Normally  about  one  4 ,5-6*0  magnitude  earthquake 
occurs  each  year  in  this  region.  About  four  other  shocks 
large  enough  to  be  felt  oc^u*  each  year.  Study  of  these 
earthquakes  has  yielded  more  reliable  regional  travel  time 
curves.  Investigation  of  ';he  region  is  hampered,  however* 
by  lack  of  seismic  stations  in  the  immediate  epicentral  re¬ 
gion.  More  recently  two  temporary  Jung-recordin*,*  trailer 
inrtslled  stations  have  been  set  up  at  epicentral  distances 
of  the  order  of  6D  km  and  less.  While  in  the  past  larger 
earthquakes  have  occurred  with  apparently  no  aftershock 
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sequence,  we  have  obst  ved  at  least  one  aftershock  series 


following  upon  two  magnitude  4  earthquakes  of  25  May  1964 


We  also  find  that  the  level  of  seismic  activity  includes 


an  average  of  one  small  earthquake  per  week.  By  placing 


untended  stations  in  boreholes  in  the  immediate  epicentrel 


region  we  feel  we  esn  for  the  first  tiir  obtain  reliable  in¬ 


formation  on  focal  depths. 


Following  a  suggestion  made  by  Stanford  Research  In¬ 


stitute,  after  an  unsuccessful  attempt  to  monitor  after¬ 


shocks  of  the  southeast  Misscur^  earthquake  of  3  March  1963, 


ARFA  sponsored  a  calibration  shot  in  the  area  to  oneck  ac¬ 


curacy  of  the  epicenter  location  and  to  test  the  calibration 


procedure.  Accuracy  of  the  epicenter  location  was  convin¬ 


cingly  checked]  the  calibration  shot,  code  named  BILLIKEN, 


was  located  within  2-3  km  of  he  shot  point  using  conven¬ 


tional  procedures  for  locating  earthquake  epicenters  in  the 


region.  The  calibration  shot  has  been  of  value  in  subse¬ 


quent  epicenter  locations  and  has  also  jeen  helpful  in  cb 


taining  some  measure  of  fecal  depth  for  the  e,  chquakes  of 


2  February  1962  and  3  March  1963 


A  joint  U.3.  Geological  Survey  -  Saint  Louis  University 


;xpiosion  refraction  survey  was  conducted  along  two  lines  in 


Missouri  in  June  1^03-  The  reduced  travel-time  plot  shows 


arrivals  to  be  better  fitted  by  a  continuously  charging  curve 


rather  than  by  a  small  number  of  straight  line  segments. 


Accordingly,  interpretation  of  the  data  has  proceeded  using 


crustal  models  characterised  by  '"ontinuous,  increasing  valoc 
ity-depth  functions.  Two  computer  programs  have  been  de- 
velopedk  one  which  calculates  travel-time  curves  and  ampli¬ 
tudes  for  a  given  velocity-depth  function;  a  second  (more 
interesting)  which  aoes  the  inverse  problem  -  it  calculates 
the  velocity  depth  function  from  a  given  travel  time  curve . 
This  work  is  unpublished  as  yet,  but  is  nearing  completion. 


2.  Abstracts 


2.1  Network  of  Seismic  Station! 


THE  RESPONSE  CHARACTERISTICS  OF  THE  LONG  PERIOD  SEISMO 
GRAPHS  OF  THE  SAINT  LOUIS  UNIVERSITY  NETWORK. 


Otto  Nuttli  and  Thomas  V.  McEvilly 
Vol.  XXXII,  Sept. -Dec.  1961,  27-36 


Earthquake  Notes 


Abstract:  Saint  Louis  University,  under  a  VELA-Uniform  con- 
trrct,  has  set  up  a  quadrilateral  network  of  permanent  seis¬ 
mograph  stations  in  the  central  United  States.  These  stations 
are  operated  with  the  cooperation  of  the  University  of  Indiana 
at  Bloomington,  the  University  of  Missouri  at  Rolla,  Kansas 
State  University  at  Manhattan,  and  Loras  College  at  Dubuque, 
Iowa,  The  latitude  and  longitude  of  the  stations  are  given 
in  Table  II  of  the  text. 


The  purpose  of  this  paper  is  to  give  a  uescription  of 
the  long -period  seismographs  operating  at  these  stations. 
Because  the  seismographs  are  intended  to  provide  data  for 
both  the  study  of  focal  mechanisms  using  S  waves  and  the 
phase  and  group  velocity  of  long-period  surface  waves,  the 
seismograph  systems  were  designed  so  as  to  give  a  nearly 
constant  magnification  for  periods  of  about  two  to  fifty 
seconds,  a  smoothly  and  slowly  changing  instrumental  phase 
shift  with  increasing  period  of  the  earth  motion,  and  sta¬ 
bility  of  operation.  The  system  adopted  uses  a  long -period 
Sprengnother  seismometer  operating  with  a  natural  period  of 
fifteen  seconds,  with  the  output  of  the  two  coils  connected 
in  series.  Shunt  and  attenuator  resistances  of  three  hun¬ 
dred  ohms  each  give  a  loose  coupling  between  the  seismometer 
and  ninety-second  galvanometer.  The  system  has  a  magnifica¬ 
tion  of  about  one  thousand  from  two  to  fifty  seconds,  and 
of  one  hundred  at  a  period  of  approximately  two  hundred 
seconds . 


Magnification  and  velocity  sensitivity  curves  are  given 
for  the  three  conqponents  of  the  long-period  Instruments  at 
each  of  the  four  stations.  In  addition,  phase  and  group 
delay  correction  curves  are  given  for  one  of  the  seismograph! 


2.2  Phase  and  Group  Velocity  Studies, 

CENTRAL  U.S.  CRUST-UPPER  MANTLE  STRUCTURE  FROM  LOVE  AND  RAY¬ 
LEIGH  WAVE  PHASE  VELOCITY  INVERSION. 

T.  V.  McEvilly,  Bull.  Seism.  Soo .  Am.,  54: 1997-2015. 

Abstract:  Fundamental  mode  Rayleigh  and  Love  wave  phase 
velocities  In  the  period  range  5  to  80  seconds  have  been 
measured  in  the  central  U.S.  Using  the  concept  of  constant 
partial  derivatives  of  phase  velocity  with  respect  to  layer 
parameters  in  a  least-squares  inversion  scheme  a  structural 
model  for  the  area  was  determined.  Several  parameters  were 
constrained  on  the  basis  of  other  available  data.  It  was 
not  possible  to  fit  satisfactorily  both  Rayleigh  and  love 
data  with  a  simple  isotropic  model.  Conclusions  require 
either  anisotropy  in  the  upper  mr.ntle  or  systematic  errors 
in  the  phase  velocity  measurements. 

EFFECT  OF  SEDIMENTARY  THICKNESS  ON  SHORT-PERIOD  RAYLEIGH 
WAVE  DISPERSION. 

T.  V.  McEvilly  and  William  Stauder,  S.J.,  Geophys.,  30, 
198-203.  — 

Abstract:  large  differences  in  group  velocities  of  short- 
perTocTRayleigh  waves  from  stripmine  blasts  for  different 
propagation  paths  in  the  Ozark  Uplift-Iliinois  Basin  area  have 
been  observed.  Good  well  control  in  the  area  makes  possible 
the  construction  of  structural  models  of  the  sediments - 
basement  system  for  these  paths.  Theoretical  group  veloc¬ 
ities  computed  for  these  models  agree  well  with  observations, 
thus  explaining  the  large  variations  in  velocities  in  terras 
of  basement-depth  differences.  This  sensitivity  of  short- 
period  surface  waves  to  sedimentary  thickness  suggests  an 
inexpensive,  single-station  technique  of  ba3in  reconnaissance 
»,..cre  coimnereial  blasting  is  available. 

DISPERSION  OF  RAYLEIGH  WAVES  IN  THE  GULF  OF  MEXICO  AND 
CARIBBEAN  SEA, 

Basil  C.  Papazachos,  Bull,  Seism.  Soc.  Am.,  54:  909-925. 

Abstract:  Group  and  phase  velocities  of  Rayleigh  waves 
along  37  paths  between  Central  America  and  the  central 
United  States  were  determined.  The  waves  examined  were  of 
period  between  15  and  60  seconds  and  crossed  the  Gulf  of 
Mexico  and  Caribbean  Sea  in  many  directions .  The  follow¬ 
ing  conclusions  about  the  structure  of  the  region  are  drawn: 

Along  the  paths  crossing  the  western  and  central  Gulf 
and  the  regions  near  the  coast  of  Texas  and  the  eastern 
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coast  of  Mexico  the  velocities  were  found  to  be  very  low. 
This  is  attributed  to  the  great  thickness  of  the  sediments 
in  these  regions. 

The-  velocities  along  paths  crossing  the  eastern  Gulf 
and  westernmost  Caribbean  were  found  to  be  higher  than  the 
typical  continental  velocities.  This  is  considered  as  an 
indication  of  a  difference  in  structure  between  the  eastern 
and  western  Gulf  of  Mexico.  Along  one  path  in  the  western 
Caribbean  Sea  between  Colombia  and  Cuba  the  velocities  were 
almost  continental.  High  velocities  found  for  part  of  the 
eastern  Caribbean  indicate  an  almost  oceanic  chaiacter  for 
this  region. 

For  each  earthquake  examined  the  initial  phases  were 
determined  with  an  error  less  than  +  tt/4.  For  a  given 
earthquake  in  all  cases  except  two  The  initial  phases  were 
found  to"  be  the  same  at  all  stations. 


2o3  Reduction  of  Ground  Motion  to  Motion  in  Incident  Wave. 
2.31  P  Waves . 

SOME  OBSERVATIONS  RELATING  TO  THE  EFFECT  OF  THE  CRUST  ON 
LONG-PERIOD  P-WAVE  MOTION. 

Otto  Nuttli,  Bull,.  Seism.  Soc.  Am.,  141-149. 

Abstract:  Long-period  P  waves  recorded  in  the  central 
United  States  at  distances  of  20  to  110  degrees  for  earth¬ 
quakes  of  magnitude  6  and  greater  Indicate  that  the  ampli¬ 
tude  of  the  ground  motion  is  ^ery  similar  at  all  stations 
in  the  area.  The  crust  is  * uiform  In  this  region,  but  the 
sedimentary  section  varies  both  in  thickness  and  in  charac¬ 
ter.  Direction  of  approach,  as  determined  by  the  ratio  of 
the  amplitudes  of  the  horizontal  components  of  motion,  does 
not  differ  significantly  from  the  great-circle  path  in  the 
majority  of  cases  studied. 

Apparent  angles  of  incidence  cannot  be  satisfactorily 
explained  by  a  half -space  model,  because  calculations  based 
on  such  a  model  require  a  larger  value  of  the  P-wave  veloc¬ 
ity  for  the  shorter  period  waves.  However,  calculations 
based  upon  Haskell’s  model  of  a  crustal  layer  over  a  half¬ 
space  give  a  satisfactory  fit  to  the  observed  data. 
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SOME  EFFECTS  OF  THE  CRUST  AND  FREE  SURFACE  ON  THE  AMPLI¬ 
TUDES  OF  BODY  WAVES. 

Otto  Nuttli,  VESIAC  Conf.  on  Variations  of  the  Earth1 s 
Crust  and  Upper  Mantle,  July  1964. 

Abstract:  Both  theoretical  and  observational  studies  have 
demonstrated  that  the  motion  recorded  by  a  seismograph,  after 
reduction  for  the  response  characteristics  of  the  instrument, 
does  not  correspond  to  the  motion  of  a  body  wave  just  before 
it  spikes  the  earth1  s  surface.  This  paper  concerns  the  ef¬ 
fects  of  the  cruLt  and  free  surface  on  data  used  in  magni¬ 
tude  and  focal  mechanism  studies,  and  the  applications  of 
such  effects  to  the  determination  of  the  crustal  structure 
at  the  recording  site. 

The  principal  conclusions  are: 

(a)  The  crustal  and  upper-mantle  structure  has  an  im¬ 
portant  effect  on  the  amplitude  and  wave  shape  of  recorded 
body-wave  motion.  However,  these  effects  depend  not  only 
on  the  crustal  structure,  but  also  upon  the  frequency  con¬ 
tent  of  the  incident  signal  and  the  epicentral  distance. 

(b)  A  consideration  of  the  effects  of  the  crust  on  the 
amplitudes  of  the  recorded  P-wave  motion  may  help  to  recon¬ 
cile  differences  observed  between  magnitudes  (m)  determined 
from  short-period  waves  and  those  from  long-period  waves. 

(c)  It  is  necessary  to  consider  the  effect  of  crustal 
structure  when  determining  the  polarization  angle  of  shear 
waves  for  epicentral  distances  less  than  about  5500 to  6000 

kiu . 


(d)  It  Is  possible,  at  least  in  theory,  to  use  observa¬ 
tions  of  recorded  body- wave  motion  to  infer  crust  and  upper- 
mantle  structure.  In  order  to  obtain  a  high  resolving  power, 
it  will  probably  be  necessary  to  work  with  frequencies 
greater  than  5  cps.  However,  if  one  is  primarily  interested 
in  total  crustal  thickness  and  average  properties,  lower- 
frequency  waves  (periods  on  the  order  of  ten  seconds)  may 
be  used.  The  latter  xrequency  range  has  the  advantage  of 
necessitating  fewer  arithmetical  calculations. 


ANGLE  OF  INCIDENCE  AND  AMPLITUDE  RATIO  OF  P  AND  PP  WAVES. 

Basil  Papazachos .  Bull.  Seism.  Soc.  Am.,  54:  105-121. 

Abstract:  Apparent  angles  of  incidence  .and  amplitudes  of  P 
and'  PP  waves  of  period  7  to  15  seconds  were  obtained  from 
seismograms  of  the  long  period  Colurabia^Sprengriether  instru 
mente  at  the  stations  Manhattan  (Kansas!,  Holla  (Missouri), 
Dubuque  (Iowa)  and  Bloomington  (Indiana).  The  calculated 


angle  of  incidence,  by  assuming  Poisson's  ratio  equal  to  £, 
agrees  with  a  velocity  of  6.4  km/sec  for  these  longitudinal 
waves  near  the  surface  of  the  earth.  T  \e  amplitude  ratios 
of  reflected  to  incident  waves  at  the  points  half  the  dis¬ 
tance  between  epicenters  and  stations  were  calculated.  CCm 
parison  of  these  ratios  with  theoretical  results  shows  that 
the  unsymmetricai  radiation  at  foci  and  errors  in  measure¬ 
ments  make  difficult,  but  not  impossible,  the  use  of  these 
ratios  for  studying  the  structure  of  the  earth.  Another 
possible  way  to  use  the  relative  amplitudes  of  P  and  PP 
waves  for  studying  regions  of  doubtful  structure  is  to  ex¬ 
amine  if  those  regions  give  reflections  not  only  at  the 
surface  of  the  earth  but  also  at  the  bottom  of  a  possible 
crustal  layer.  Atlantic  reflections  indicate  that  these 
longitudinal  waves  are  not  affected  by  any  surface  crustal 
layer  under  the  Atlantic  Ocean. 


AN  APPLICATION  OP  THE  HASKELL-THOMSON  MATRIX  METHOD  TO  THE 
SYNTHESIS  OF  THE  SURFACE  MOTION  DUE  TO  DILATATIONAL  WAVES. 


W.  J.  Hannon,  Bull.  Seism.  Soc.  Am 


Abstract  The  surface  motion  arising  from  a  plane  dilatation 
ai  wave  striking  the  base  of  a  layered  system  is  obtained  by 
Fourier  synthesis  using  the  frequency  dependent  transmission 
coefficients  which  may  be  computed  by  the  Haskell -Thoms on 
matrix  method.  Theoretical  seismograms  il1  strating  the 
effect-  of  thin,  low- velocity  surface  layers  on  the  surface 
motion  are  constructed  from  the  transmission  coefficients 
of  four  crustal  models. 


2.32  S  Waves 


ON  THE  DETERMINATION  OF  THE  POLARIZATION  ANGLE  OF  THE  S  WAVE 


Otto  Nuttli  and  John  D.  Whitmore,  Bull.  Seism*  Soc.  Am 
52:  95-107. 


Abstract:  This  study  is  concerned  with  determining  the  mini 


mum  epicentral  distance  for  which  it  is  permissible  to  ob¬ 
tain  the  value  of  the  polarization  angle  of  the  S  wave  by 
measuring  the  angle  between  the  great  circle  path  at  the 
station  and  the  direction  of  the  horizontal  component  of 
the  S  wave  particle  motion  obtained  from  the  seismograms. 
This  critical  distance  can  be  determined  by  Ine  fact  that 
at  smaller  distances  the  particle  motion  of  the  earth's 
surface  due  to  the  incidence  of  S  will  be  nonlinear  (the 
SH  and  the  horizontal  and  vertical  components  of  SV  will  be 
out  of  phase  with  respect  to  one  another)  while  at  larger 


distances  the  particle  motion  v/ill  be  linear.  An  analysis 
of  the  S  motion  recorded  by  the  Galitzin-ltfilip  seismographs 
at  Florissant  indicates  that  the  critical  distance  is  42 
degrees.  The  periods  of  these  S  waves  are  of  the  order  of 
10  second.  The  analysis  also  shows  that-  the  effective  P 
wave  velocity  of  teleseismic  waves  at  the  earth's  free  sur¬ 
face  is  7,74  kis/sec,  and  the  effective  value  of  Poisson's 
ratio  and  the  effective  S  wave  velocity  at  the  earth's  sur¬ 
face  are  0.25  and  4.46  km/sec,  respectively.  By  effective 
values  are  meant  the  values  of  the  velocities  and  Poisson's 
ratio  that  govern  the  angle  of  incidence  of  the  waves  at  the 
earth's  surface. 


THE  DETERMINATION  OF  S-WAVE  POLARIZATION  ANGLES  FOR  AN 
EARTH  MODEL  WITH  CRUSTAL  LAYERING. 

Otto  Nuttll,  Bull.  Seism,  Soc.  Am.,  £4:  1429-1440„ 

Abstract:  This  paper  presents  a  method  for  determining  the 
polarization  angle  of  S  waves  which  takes  account  of  the 
crustal  layering  at  the  seismograph  station.  Charts  are 
given  for  four  crustal  models,  corresponding  to  normal  con¬ 
tinental,  thick,  and  thin  crusts,  which  enable  one  to  ob¬ 
tain  the  polarization  angle  at  the  top  of  the  mantle  beneath 
the  seismograph  station. 

The  equations  to  be  used  for  obtaining  the  polariza¬ 
tion  angle  are 

tan  £  «  Fh  tan  y 

tan  £  =  Fv  tan  £' 

where  £  is  the  polarization  angle,  y  is  the  angle  between  th 
horizontal  component  e**  the  S-wave  surface  motion  and  the 
great-  circle  path  at  the  station,  &  is  the  angle  oetween  the 
vertical  axis  and  the  component  of  the  S-wave  surface  motion 
in  the  plane  transverse  to  the  great  circle  path,  and  F^  and 
Fv  are  functions  cf  wave  period,  eplcentral  distance,  and 
the  crustal  structure  at  the  seismograph  station.  Values  of 
Fh  and  Fv  for  four  crustal  models  are  given  in  the  paper. 


A  METHOD  FOR  DETERMINING  THE  TIME  OF  ONSET  OF  S  WAVES  OF 
SHALLOW  EARTHQUAKES  AT  INTERMEDIATE  DISTANCES. 

Otto  Nuttll  and  Ahou-Bakr  Ibrahim,  Paper  presented  at  the 
Annual  Meeting  of  the  Seismological  Society  of  America, 
Seattle,  Washington,  March  28-30,  1964. 

Abstract:  Observational  data  concerning  the  travel  times 


of  S  waves  from  shallow  earthquakes  at  distances  of  about 
2000  to  4000  kilometers  are  very  limited,  owing  to  the 
difficulty  in  identifying  the  beginning  of  the  S-wave  mo¬ 
tion.  Both  interfering  longitudinal  waves  and  phase 
shifts  in  S  due  to  crustal  layering  and  the  earth's  free 
surface  obscure  the  onset  of  the  S  arrival.  However,  a 
comparison  of  observed  with  theoretical  particle  motion 
diagrams  enables  one  to  determine  the  beginning  of  S  in 
cases  where  it  is  difficult  to  do  so  from  a  visual  inspec 
tior.  of  the  seismograms.  Examples  of  the  application  of 
this  method,  including  theoretical  particle  motion  dia¬ 
grams  for  an  earth  model  with  crustal  layering,  are  pre¬ 
sented  . 


2.4  Focal  Mechanism  Studies 


SOME  P  AND  S  WIVE  STUDIES  OF  THE  MECHANISM  OF  EARTHQUAKES 
AND  SMALL  CHEMICAL  EXPLOSIONS. 


William  St&uder,  S.J.,  and  Car]  Kisslinger,  VESIAC  Spec 
Rpt,  Proceedings  of  the  Colloquium  on  Detection  of  Under 
ground  Nuclear  Explosions,  August  1962,  pp. 125-140. 


Abstract 


Both  the  earthquake  and  explosion  studies  point 
.  for  further  research  into  the  nature  of  seismic 
source  ffiechanisms  as  such.  They  show  that  careful  observa¬ 
tions  of  S  waves  are  a  useful  tool  in  this  type  of  investi¬ 
gation-  and  that  the  S  waves  are  diagnostic  of  the  pliysical 
processes  taking  place  at  the  source. 


As  yet,  no  difference  in  the  S  waves  from  the  two  types 
of  sources  has  been  discovered  that  will  indicate  which 
source  is  operative.  The  presence  or  absence  of  S  waves 
is  certainly  not  a  criterion.  The  ratio  of  S  to  P  or  of 
Rayleigh  to  Love  may  be  indicative,  but  much  research  re¬ 
mains  to  be  done  in  this  area.  Moreover,  such  studies 
must  not  ignore  the  effect  of  the  medium  surrounding  the 
explosion  on  the  amplitude  of  S  waves  relative  to  P. 

Unfortunately,  data  from  nuclear  explosions  are  avail¬ 
able  for  a  limited  variety  of  media.  We  suggest  that 
future  underground  nuclear  explosions  be  given  the  most 
complete  azlrautha]  coverage  practical  and  that  this  cover¬ 
age  be  with  three  consonant  instruments.  In  addition, 
careful  pre-  and  postshot  surveys  cf  cracks  and  the  rela¬ 
tion  of  these  tc  regional  joint  patterns  should  be  included 


S-WAVE  STUDIES  OF  EARTHQUAKES  OF  THE  NORTH  PACIFIC,  Pert  I 
KAMCHATKA 


William  Stauder,  S.J.,  Bull.  Seism*  So 
527-550. 


Abstract;  The  polarization  of  S  waves  at  stations  di&- 
triBufed  a zimut hally  about  the  source  is  examined  for  eaeh 
of  twenty-three  Kamchatka  earthquakes  of  195C-1960,  In 
nineteen  of  these  earthquakes  the  P  and  S  wave  data  are  in 
agreement  with  a  double  couple  source  as  the  point  model 
of  the  focal  mechanism.  The  S  waves  indicate  a  uniform 
mechanism  which  repeats  itself  from  earthquake  to  earth¬ 
quake  and  from  which  It  may  be  inferred  that  the  axes  of 
greatest  and  least  stress  at  the  foci  tend  to  lie  in  a 
vertical  plane  normal  to  the  trend  of  the  Kamchatka-Kuriles 
arc.  The  axis  of  least  stress  usually  plunges  almost  ver¬ 
tically-  under  the  continent,  but  may  also  plunge  less 
steeply,  at  angles  as  low  as  4?  degrees*  At  least  two 
earthquakes  may  be  represented  by  a  single  couple  source. 


S- WAITE  STUDIES  OF  EARTHQUAKES  OF  TK  '  NORTH  PACIFIC.  FART  II: 
ALEUTIAN  ISLANDS. 


William  Stauder,  S.J.,  and  Agustin  Udia&,  S.J.,  Bull 
Seism.  Soc.  Am.,  53:  59-77. 


Abstract:  The  polarization  of  the  S  wave  at  stations  dis¬ 
tributed  azinruthally  about  the  source  is  examined  for  each 
of  twenty-five  Aleutian  Island  earthquakes.  A  combination 
of  data  from  the  first  motion  of  P  and  from  the  polariza¬ 
tion  of  S  is  then  used  to  study  the  focal  mechanisms  of  the 
earthquakes.  This  combination  of  P  and  S  wave  data  is  found 
to  make  possible  a  good  determination  of  the  focal  mechan¬ 
ism  in  eases  where  data  from  the  first  motion  of  P  alone 
do  not  suffice. 


The  earthquakes  are  divided  into  three  groups  according 
to  three  basic  patterns  of  S  wave  polarization.  The  first 
group  (fourteen  earthquakes)  corresponds  to  a  double  couple 
The  second  group  (five  earthquakes ) and  the  third  group  (six 
earthquakes)  are  conformable  to  conjugate  shears  and  may 
therefore  be  explained  by  single  couple  sources  of  opposite 
moment,  respectively.  It  is  shown  that  a  uniform  principal 
stress  system  predominates  in  the  region  and  that  the  axis 
of  greatest  compressive  stress  is  normal  to  the  trend  of 
the  island  arc. 
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A  LEAST  SQUARES  METHOD  FOR  EARTHQUAKE  MECHANISM  DETERMINA¬ 
TION  USING  S-WAVE  DATA. 


Ague  tin  Udias,  S.J.,  Ball.  Seism.  Soc.  Am.,  54,  2037-2047. 


Abstract;  In  this  paper  a  numerical  approach  to  the  deter¬ 
mination  of  focal  mechanisms  based  on  the  observation  of 
the  polarization  of  the  S  wave  at  N  stations  is  presented. 
Least -suqare  methods  are  developed  for  the  determination 
of  the  orientation  of  the  3ingle  and  double  couple  sources, 
The  methods  allow  a  statistical  evaluation  of  the  data  and 
of  the  accuracy  of  the  solutions. 


APPLICATION  OF  NUMERICAL  METHOD  FOR  S-WAVE  FOCAL  MECHAN¬ 
ISM  DETERMINATIONS  TO  EARTHQUAKES  OF  KAMCHATKA -KURILE 
ISLANDS  RECION. 


Agustin  Udlas,  S.J.,  and  William  Stauder,  S.J.,  Bull. 
Seism.  Soc.  Am.,  54,  2049-2065, 


Abstract:  The  least  squares  method  for  the  determination  of 
Tocai mechanism  using  S-wave  data  developed  by  Udias  (1964) 
has  been  applied  to  34  earthquakes  of  the  Karachatka-Kurile 
Islands  Regions,  For  most  of  the  earthquakes  studied  the 
source  mechanism  was  found  to  conform  best  to  a  double 
couple  model.  The  accuracy  and  stability  of  the  solutions 
show  a  steeply  dipping  axis  of  tension  and  an  almost  hori¬ 
zontal  axis  of  pressure.  The  orientation  of  the  pressure 
axes  is  normal  to  the  regional  trend.  The  direction  of 
the  stresses  suggest  a  thrust  faulting  with  predominantly 
dip  slip  motions  such  that  the  continental  block  moves 
over  the  oceanic. 


GEOMETRICAL  PROPERTIES  OF  GROUPS  OF  FAULT  PLANE  SOLUTIONS 


Allan  Cox  and  William  Stauder,  S.J.,  Bull.  Seism.  Soc.  Am. 
£4:  87-104.  - - 


Abstract:  A  formal  method  of  analysis  is  proposed  whereby 
groups  of  first  motion  solutions  may  be  examined  for  evi¬ 
dence  of  regional  patterns  of  symmetry,  abstracting  from 
any  particular  model  of  the  focus  or  from  assumptions  about 
regional  tectonics ,  Axial,  planar,  and  conical  symmetry 
are  discussed  with  special  consideration  being  given  to  the 
constraints  placed  on  the  special  distribution  of  axes  by 
the  geometry  of  sets  of  orthogonal,  axes.  Examples  from 
the  literature  are  given  of  axial  and  planar  symmetry,  and 
it  Is  shown  that  in  some  instances  attention  to  the  geome¬ 
try  of  the  axis  system  suggests  modifications  in  the  inter¬ 
pretation  of  regional  r>atterns.  Analytical  aids  are  ilso 
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given  to  associate  the  use  of  eigenvalues  with  the  graphi¬ 
cal  netnods  in  discerning  the  type  of  symmetry  and  indi¬ 
cating  the  degree  of  scatter,  and  a  method  is  described  for 
calculating  statistical  confidence  limits  for  the  case  of 
axial  symmetry. 


A  COMPARISON  OP  MULTIPLE  SOLUTIONS  OF  FOCAL  MECHANISMS. 

William  Stauder,  S.J.,  Bull.  Seism,  Soo.  Am.,  54,  927-937. 

Abstract:  The  occurrence  of  focal  mechanism  solutions  of 
TKe  same  earthquake  by  different  authors  is  examined.  Of 
42  instances  of  such  multiple  solutions  for  earthquakes 
which  occurred  around  the  borders  of  the  Pacific  Ocean, 
it  is  found  that  of  84  possible  pairs  of  multiple  solu¬ 
tions  agreement  is  excellent  in  8  cases,  good  in  13,  fair 
in  19,  poor  in  19,  and  no  agreement  in  25.  The  implica¬ 
tion  for  focal  mechanism  studies  rs  illustrated  by  examin¬ 
ing  the  focal  mechanism  solutions  of  deep-focus  earth¬ 
quakes  under  the  Sea  of  Japan. 


THE  FOCAL  MECHANISM  OF  TWO  PERUVIAN  EARTHQUAKES, 

Daniel  E,  Cchoa  and  William  Stauaer,  3.J.,  Pure  &  Applied 
geophysics,  58  (1964 /II):  114-122.  - 

Abstract:  The  focal  mechanisms  of  two  Peruvian  earth¬ 
quakes,  January  15,  1958,  and  January  13,  i960,  are  deter¬ 
mined  using  data  from  both  the  first  motion  of  P  and  the 
polarization  of  S,  The  fault  motions  correspond  to  strike 
slip  motion  on  neighboring  faults  corresponding  in  strike 
to  geophysical  evidences  of  regional  faulting.  The  motion 
is  left  lateral  in  one  of  the  shocks,  right  lateral  in  the 
other . 


THE  PRESENT  STATUS  OF  FAULT  PLANE  SOL^IONS  FROM  P  AND  S 
WAVES  AND  THE  RELATION  OF  THESE  SOLUTIONS  TO  TECTCNIC 
STRESS  FIELD. 

William  Stauder,  S. J.  Proceedings  VESIAC  Conference  on 
Focal  Mechanisms,  LaJolla,  March  1965,  (in  press)  Univer¬ 
sity  jf  Michigan. 

Abstract^  The  P  wave  first  motion  method  of  focal  raschan- 
Isni  de Fermi na t ion  is  used  to  infer  the  symmetry  of  a  seis¬ 
mic  source.  In  general,  three  classes  of  solutions  may  be 
distinguished:  those  in  which  two  P  nodal  planes  are 
clearly  determined,  those  in  which  only  one  nodal  plane 


is  clearly  determined,  those  in  which  neither  nodal  plane 
is  distinguished,  i.e.,  all  available  first  motion  readings 
are  of  the  same  sign.  Even  in  the  latter  case,  while  the 
plane  of  faulting  or  dir*  tion  of  motion  cannot  be  discerned 
considerable  information  about  the  stress  field  at  the 
source  of  the  shock  is  indicated;  for  instance,  if  all  first 
motion  readings  are  congressional,  the  axis  of  least  prin- 
cipal  stress  is  steeply  inclined. 

Inability  to  aeterrine  both  nodal  planes  in  all  cases 
results  from  a  limitation  inherent  to  the  method,  that  is, 
rays  leaving  the  source  and  arriving  at  stations  at  Pn 
distances  or  greater  all  leave  the  focus  downwa  d  through 
a  narrow  cone  which  represents  3/4  or  less  of  the  area  of 
the  focal  sphere.  A  related  limitation  arises  from  the  geo¬ 
graphic  distribution  of  ...  ;ismic  stations.  These  and  other 
limitations  emphasize  the  desirability  of  using  S  wave 
data  and  r  urface  wave  data  in  focal  mechanism  determina¬ 
tion.  Methods  of  i  vave  analysis  and  surface  wave  analysis 
hav5,  in  face,  been  developed  and  have  been  used  success¬ 
fully  in  focal  mechanism  studies.  These  methods  have 
greatly  enhanced  the  reliability  of  fault  plane  solutions. 

In  examining  groups  of  fault  plane  solutions,  consider¬ 
able  evidence  of  bias  in  the  solutic.is,  of  unreliability 
of  individual  solutions,  and  of  systematic  differences  in 
solutions  by  different  authors  Indicate  that  caution  should 
be  exercised  in  undertaking  statistical  studies  of  all  the 
fault- plane  solutions  for  a  given  region.  Nevertheless, 
where  good  solutions  are  available,  consistent  relations 
between  solutions  are  discerned.  Most  notably,  around  the 
borders  of  the  North  Pacific  the  axes  of  greatest  and  I  Bast 
stress  seem  to  lie  in  a  steeply  dipping  plane. normal  to 
the  trend  of  the  geographic  features.  The  axis  cf  coraore^ 
sive  stress  tends  to  be  inclined  hy  an  angle  of  10°  to  30° 
to  the  horizontal. 


THE  DIRECTION  OP  CRUSTAL  MOVEMENTS  INDICATED  BY  EARTH 
QUAKE  DATA 


William  Stauder,  S.J.  Paper  delivered  at  the  Research 
Symposium,  Annual  Vfr  3 ting  of  the  AAPC-SEPM,  New  Orleans 
April  26-29,  1965. 


Absto  ct;  Two  types  of  seismic  data  give  evidence  of  the 
direction  of  crustal  movements:  the  spatial  ana  temporal 
sequence  of  earthquake  Lypocenters  and  the  focal  mechan¬ 
ism  determinations  of  individual  earthquakes  or  of  groups 
of  earthquakes.  Examples  illustrate  the  power  of  the 
methods  of  focal  mechanism  determination  irom  P  and  S  waves 


as  also  the  limitations  inherent  to  the  methods.  .Further 
examination  of  groups  of  published  fault  plane  solutions 
shows  evidence  that  many  of  these  solutions  are  not  relia 
ble.  A  critical  examination  of  the  re?,iability  of  in¬ 
dividual  fault  plane  solutions  before  these  are  used  as 
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Conclusion  of  tectonic  significance  from  the  data  of 
earthquakes  include  the  following.  1  rst,  the  predom¬ 
inance  of  strike  slip  faulting,  once  reported  as  a  con¬ 
clusion  of  focal  mechanism  studies,  is  no  longer  apparenc 
from  the  data  of  fault  plane  solutions.  Second,  around 
the  borders  of  the  northern  and  western  Pacific  Ocean  the 
predominant  horizontal  stress  is  a  compression  normal  to 
the  trtr.d  cf  the  tectonic  features.  Third,  there  is  a 
radical  difference  between  the  earthquake  zones  of  the 
circum-Paclfic  or  Alpide  be3.t,  and  these  of  the  mid-oceanic 
rifts;  the  latte**  represent  very  narrow  zones  of  seismic 
activity  in  which  the  seismic  evidence,  as  yet  incomplete, 
indicates  that  the  rift  is  indeed  a  tensional  feature, 
fourth,  that  along  the  world  rift  system  concentrations  of 
epicenters  are  associated  with  offsets  of  the  oceanic 
ridges.  These  offsets,  in  turn,  correspond  to  the  inter¬ 
section  of  the  ridge  by  east-west  trending  fracture  zones. 
Along  the  fracture  zones  themselves  the  current  seismic 
activity  is  limiteo  to  the  region  of  offset.  These  rela¬ 
tions  have  been  especially  invesuigated  by  Sykes.  Fifth, 
Sykes  has  also  shown  that  in  the  island  arc  structure  of 
the  Kermadec -Tonga  Islands,  earthquake  hypocenters  lie 
along  a  pi a  f  dipping  Inward  below  the  island  arc,  and 
that  at  all  levels  the  hypocenters  reflect  the  surface 
details  of  the  arc.  Evidence  from  the  focal  mechanism 
solutions  of  earthquakes  occurring  along  the  arc  suggest 
that  the  planes  of  faulting  fellow  this  same  pattern. 


Focal  Depth 


THE  TRANSFER  FUNCTION  FOR  P-WAVES  FOR  A  SYSTEM  CONSISTING 
OF  A  POINT  SOURCE  If  A  LAYERED  MEDIUM. 


Karl  Fueh3 ,  Scientific  Report  No.  12  prepared  under 
Contract  AF  19(604} -7399  of  Air  Force  Cambridge  Research 
Laboratories  for  ARTA  Project  VELA-Uniform,  AFORL-65-526 


Abstract:  This  paper  investigates  In  what  manner  the  spec¬ 
trum  ofbody  waves  radiating  from  point  sources  in  a 
multilayered  medium  ow*  a  homogeneous  half-space  is  dif¬ 
ferent  from  the  spectrum  of  body  waves  frGm  the  same  source 
in.  an  infinite  medium.  The  effect  of  the  system  consisting 
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of  a  point  source  in  a  layered  crust  on  the  spectrum  of  p- 
waves  observed  at  large  distances  in  the  half  space  is 
studied.  Analytical  expressions  for  Uie  transfer  function 
of  this  system  are  derived  for  three  types  of  point  sources: 
an  explosive  source,  a  single  couple,  and  a  double  couple 
of  arbitrary  orientation  within  the  crust. 

Preliminary  numerical  computations  for  the  explosive 
source  at  various  depths  in  a  realistic  model  of  the  earth’s 
crust  study  the  effect  of:  a)  the  angle  of  incidence  Into 
the  homogeneous  half-space,  b)  the  source  depth,  c)  minor 
variations  of  the  crustal  model.  In  the  case  of  an  explo¬ 
sive  source  the  most  influential  parameter  of  the  transfer 
function  is  the  source  depth.  Tn  shallow  explosions  t be¬ 
low  frequency  part  of  the  spectrum  of  body  waves  is  compara¬ 
tively  rejected. 


2,5  Inversion  of  Body  Wave  Spectra, 


SOME 


EFFECTS 


OF  A  LAYERED  SYSTEM  ON 


DILATATION!. 


Willard  J.  Hannon.  Technical  Report  3,  July  1964*  prepared 
for  Air  Force  Cambridge  Research  Laboratories,  OAR,  under 
ARPA  Project  VELA -Uni form.  AFCRL  Rpt.  No.  64-614. 


Abstract:  In  this  report,  the  effect  of  the  crustal  model 

on~IxiG^varIatior  of  the  surface  motion  with  the  angle  of 
incidence  and  the  frequency  is  examined  i or  several  crustal 
models.  The  study  was  carried  out  by  programming  the  prob¬ 
lem  for  the  IBM  1620  and  7072  computer  systems  using  the 
matrix  formulation  originally  suggested  by  Thomson  and  per¬ 
fected  by  Haskell  and  Dorman.  From  these  programs,  the  ratios 
of  the  displacements  at  the  free  surface  to  the  total  ampli¬ 
tude  at  depth  were  computed  for  several  crustal  models  in 
ranges  of  frequency  and  angle  of  Incidence  of  interest  in 
seismology.  These  ratios  are,  in  effect,  transmission  co¬ 
efficients. 


Six  crustal  models  having  such  features  as  thin  low- 
velocity  surface  layers,  low~veiocj.ty  layers  at  depth,  and 
relatively  thick  and  relatively  thin  total  thicknesses  were 
considered.  For  each  model,  the  transmission  coefficients 
were  computed  for  frequencies  ranging  from  .02  ops  to  10.0 
cps  in  steps  of  .02  cps,  and  for  angles  of  Incidence  ranging 
from  21  degrees  to  53  degrees  in  steps  of  4  degrees.  Has¬ 
kell’s  model  was  Included  in  these  calculations  in  order  to 
obtain  a  check  on  the  calculations.  A  further  check  was 
obtained  by  using  the  transmission  coefficients  to  synthesize 
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the  surface  motion  due  to  an  incident  wave  of  the  form 
(1/  Tf  )  (sin  SO  tft)/t  and  comparing  these  values  with  'hose 
predicted  by  ray  theory*  The  agreement  was  very  good. 

As  a  result  of  these  calculations,  the  importance  of 
the  frequency  dependent  character  of  the  crustal  effect  has 
been  further .emphasized.  It  has  been  shown  that  a  thin 
low-veiocity  surface  layer  causes  a  large  variation  in  the 
transmission  coefficients,  while  a  low-velocity  layer  at 
depth  has  little  effect,  especially  at  low  frequencies. 
Further,  it  has  been  shown  that  the  total  crustal  thickness 
is  one  of  the  most  significant  factors  in  determining  the 
variation  with  frequency.  The  transmission  coefficients  of 
relatively  thick  crustal  models  have  a  much  more  rapid  fre¬ 
quency  variation  then  those  of  the  relatively  thin  cru3t. 

At  frequencies  less  than  .2  cps,  Jt  is  very  difficult  to 
distinguish  between  crustal  models  of  very  nearly  the  same 
thickness. by  means  of  the  frequency  variation  of  the  trans¬ 
mission  coefficients  unless  the  internal  structure  is  very 
different . 


THE  DETERMINATION  OF  CRUSTAL  THICKNESS  FROM  THE  SPECTRUM 
OF  P  WAVES . 

Luis  M*  Fernandcx.  S.J,  Scientific  Report  No.  13,  Sept. 
1965*  prepared  for  hir  Force  Cambridge  Research  Laboratories, 
OAR,  under  Contract  AF  19(604)-7399,  ARP  A  Project  VELA- 
Uniforra. 

Abstract:  When  the  layers  of  the  earth1 s  cru3t  are  perturbed 
TSy^corapresaional  seismic  energy  they  vibrate  in  a  frequency 
characteristic  pattern.  This  pattern  depends  not  only  on 
the  frequency  content  of  the  exciting  seisndc  energy  tut 
also  on  the  elastic  parameters  and  thicknesses  of  the  layers. 
This  dependence  may  be  on  the  basis  of  a  method  to  investi¬ 
gate  the  crust  of  the  earth. 

In  order  to  obtain  Information  independent  of  the 
source  of  energy  the  spectrum  of  the  vertical  e;  ponent  of 
motion  should  be  divided  by  the  spectrum  of  the  horizontal 
component.  This  ratio  represents  the  tangent  of  the  appar¬ 
ent  angle  of  emergence  as  a  function  of  frequency  and  as 
such  depends  only  on  the  angle  of  incidence  of  the  ray  and 
the  system  of  layers  below  the  station  of  obsf  nmtion. 

The  parameters  of  the  crust  may  be  obtained  by  compari¬ 
son  of  theoretical  and  observed  spectra.  To  facilitate  tbl3 
comparison  a  set  of  master  curves  was  calculated  using  the 
matrices  development  of  Haskell.  Calculations  of  these 
curves  are  in  terras  of  a  dimensionless  parameter  instead  of 
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frequency.  This  parameter  allows  the  grouping  of  the  curves 
corresponding  to  different  crustal  models  into  families  of 
curves.  A  set  of  master  curves  of  the  apparent  angle  of  * 

emergence  for  one-layer  models  and  for  different  angles  of 
incidence  and  contrasts  of  velocities  between  the  crust  and 
the  mantle  is  presented.  This  set  is  complete  in  the  sense 
that  all  the  one-layer  models  may  be  interpolated.  A  sec¬ 
ond  set  for  some  combinations  of  two-layer  models  is  also 
presented.  The  characteristics  of  these  curves  are  dis-~ 
cussed  from  the  point  of  view  of  periodicity  and  amplitude 
to  investigate  the  Influence  of  the  layer  parameters  on 
their  form.  Constructive  interference  considerations  and 
Fourier  analysis  of  a  pulse  multiply  reflected  within  the 
crust  reveal  that  the  amplitude  of  the  curve  fluctuations 
depends  on  the  velocity  contrast  at  the  Interfaces  of  the 
crust  and  of  the  crust-mantle.  The  periodicity  of  the  fluc¬ 
tuations  depends  on  the  time  lags  between  the  first  arrival 
of  the  direct  P  wave  and  the  secondary  arrivals  of  the  reflec¬ 
tions  and  seismic  conversions.  Long-period  fluctuations 
correspond  to  short  time  lags  and  short-period  to  long  time 
lags. 

Observations  of  the  spectrum  of  the  apparent  angle  of 
emergence  were  obtained  by  dividing  the  smoothed  spectra 
of  the  vertical  and  horizontal  component  seismograms.  In 
order  to  avoid  the  influence  of  reflections  at  the  crust 
near  the  source  or  reflections  from  the  core  of  the  earth 
the  earthquakes  selected  were  of  intermediate  and  large  fo¬ 
cal  depth. 

Application  of  the  method  to  the  long-period  seismo¬ 
grams  of  the  Saint  Louis  University  Network  of  stations 
gives  an  average  P  velocity  in  the  crust  of  6.6  Kn/sec. 
and  a  total  thickness  of  the  crust  of  42  Km,  for  the  cen¬ 
tral  part  of  the  United  States  under  these  stations. 


2.7  Effect  of  Low  Velocity  Channel  on  Travel  Time.:  of 
Body  Waves 

TRAVEL-TIME  CURVES  FOR  A  LOW-VELOCITY  CHANNEL  IN  THE 
UPPER  MANTLE. 


John  Dowling  and  Otto  Nuttil,  Bull.  Seism.  Soc.  Am..  54. 
1981-1996. 

Velocities  within  the  earth  can  be  determined  from  body 
wave  time-distance  (T-D)  data  by  the  Herglotz-WIechert 
method  provided  the  velocity  does  not  deer  ase  too  rapid¬ 
ly  with  depth.  Until  the  present  time,  the  properties  of 
T-D  curves  for  rapid  decreases  of  velocity  with  depth  have 
been  considered  only  qualitatively. 
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This  paper  presents  a  technique  for  calculating  a  T-D 
curve  for  any  velocity  distribution,  including  continuous 
and  discontinuous  increases  and  decreases  of  velocity  with 
depth.  Some  properties  of  T-D  curves  are  quantitatively 
studied  by  systematically  yarying  the  characteristics  of 
a  single  model  and  noting  the  corresponding  variations  in 
the  calculated  T-D  curves.  Prom  this  it  is  concluded  that 
a  significant  1 c w -velocity  channel  may  not  be  evidenced  by 
a  shadow  zone  but  rather  by  an  overlapping  of  two  distinct 
branches  of  the  T-D  curve.  It  is  further  concluded  that 
the  presence  of  a  shadow  zone  implies  a  very  gentle  veloc¬ 
ity  gradient  below  the  low-velocity  channel . 

By  fitting  a  calculated  T-D  curve  to  observed  data  one 
can  determine  velocity  as  a  function  of  depth  even  when  the 
velocity  decreases  rapidly  with  depth,  when  a  low-velocity 
channel  exists.  Observed  T-D  data  for  two  underground 
nuclear  explosions  (GNOME  and  BIT  BY)  measured  in  four  dif¬ 
ferent  azimuths  were  fitted  with  T-D  curves  calculated  for 
assumed  velocity  distributions .  It  is  concluded  that  these 
data  can  be  satisfied  by  a  low-velocity  channel  for  P  waves 
in  the  upper  mantle.  The  character  of  th4s  channel  (depth, 
thickness  and  velocity)  was  determined  in  each  azimuth.  The 
depth  to  its  top  was  shallow  (70  +  km)  in  the  western  U.S. 
and  deep  (125  +  km)  in  the  eastern  U.S.  The  velocity 
gradient  below  the  channel  is  sharp  enough  to  produce  no 
prominent  shadow  zones.  There  are  significant  lateral 
changes  in  upper  mantle  velocities  in  the  western  U.S. 


2.8  Local  Crus.tal  Structure,  Etc. 


Pn  VELOCITY  AND  OTHER  SEISMIC  STUDIES  FROM  THE  DATA  OF 
RECENT  SOUTHEAST  MISSOURI  EARTHQUAKES. 

lilliam  Stauder,  S.J.,  and  Gilbert  Bollinger.  Bull,  Seism. 
Soc.  Am.,  53:  661-679 

Abstract:  Five  earthquakes  occurred  in  the  southeastern  Mis¬ 
souri  seismic  region  during  the  first  seven  months  of  1962. 
The  largest  of  these  took  place  on  2  February  1962,  and  was 
located  near  New  Madrid.  The  felt  area  of  this  shock 
covered  35* OCX)  square  miles;  the  maximum  intensity  was 
Intensity  V  in  the  Modified  Mercalli  scale.  The  magnitude 
is  estimated  to  have  been  to  4  3/4. 

Pn  arrival  times  from  thirty-five  stations  in  the  dis¬ 
tance  range  20«  km  to  1000  km  fit  a  composite  travel  time 
curve  for  all  five  earthquakes  given  by 
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t  -  5.56  +  0.21  sec  +  4/(8. 23  +  0.02  km/sec). 

Phis  contrasts  with  interval  veloc-’tipe?  of*  r  ?  / 

8.5  km/sec  obtained  by  Herrin  and  Taggart  SSC  ?° 

from  the  GNOME  data  q  Qr.j  na  a,si?^art  for  this  region 

aee  also  studied?  n  P  md  3  crustal  velocities 
Pulse-like  Rayleigh  and  Love  wsvps  W0r>0  - 

is  5£FjS  k='  sSsT 

s‘^'£HSL» 

THE  SE  MISSOURI  EARTHQUAKE  OP  3  MARCH  1963 . 

N^am3f(Sept:-&c:'l|63)fi?4!el  Goodwln-  Sarthguake 

itfr^flF  a  LIs?Vel  tiroes  for*  the  earthquake  of  3  March 
relation  lea8t~sc^ai*es  travel  time  curve  given  by  the 

t  =  (6.29  +  0.29)  sec  +  A/(8.24  +  0.03  km/sec) 

Thi3peq!!Sl0^ap?lies  ^  the  distance  range  200  to  1000  km 
The  Fn  velocity  from  refraction  surveys  and  from  the  rtt™ 
experiment  are  slightly  lower.  *  Q  Ir°m  the  BILljIKEN 

lndlcate  a  maximum  intensity  of  VII 
Radius  of  the  felt  area  was  at  least  200  miles 

range  from  6  to  '+•  A  preliminary  dilc^s^ 
sion  v;ill  be  given  of  the  depth  of  focus*  focal  rre  hanism 
and  group  and  phase  velocities  of  the  sedimentary  bSh 
of  the  Love  and  Rayleigh  wave  curves. 

THE  BILLIKIN  CALIBRATION  SHOT  IN  SOUTHEAST  MISSOURI. 

John  °?Kll71S  and  Wayne  H.  Jackson. 
tl1  UT  u  RePO**t  4  prepared  for  Air  Force  Cambridge  Re- 

0AR*  under  Contract  A F  19(604^-7399 
for  ARPA  Project  VELA-Uniform.  AFCRL  64-659.  '  ^ 


^g|rgg|;.„0n  fu"e„28'  1563.  a  20,000  lb.  high-explosive 
calibration-shot,  code  named  BILLIKEN,  was  detonated  near 

h.  epicenter  01  the  magnitude  5  Southeast  Missouri  earth¬ 
quake  of  Ha  oh  3,  1963.  The  experiment  was  suggestedby 

set^n1"^  ?LS!af0r?  R?3esJ0h  Inatltl;te  after  field  units 
seo  up  in  the  epicentral  area  failed  to  record  aftershocks. 


36 


Ibe  shot  was  sponsored  by  the  Advanced  Research  Prelects 

of6the  undSC°f£n?oe?  J?1oltly  by  the  Studies ^Branch 

oi  the  U.  S,  Geological  Survey  and  Saint  Louie  University. 

Using  arrival  times  at  stations  comparable  to  those 

epleenter°oftBlLl,T^t;Mar°hi  ^  “^,‘the  instru=l 
epicenter  of  BILLIKEN  wan  located  within  2  to  4  kilometers 

ma  uie  snov  point,  Tnig  confirms  the  accuracy  of  eDicenter 

TheaI'J°"SKin  \US  area  fro“  tlie  *»ta  of  near-reglonals  tat  ions 
2'!??*  alf°  Permitted  "calibration"  of  the  area,  maKing 
S  t  locations  using  time  differences  of  S- 

Of  tho  fL?1.  Of  stations.  To  the  ea3t,  north  and  west 
thnn* 5  ®hot-P°lnt  travel  times  from  BILLIKEN  agree  with 
££???, °f  Previous  earthquake  studies  and  of  a  current  re¬ 
tain^  a SUrVeyJ i0  thf  so,lt/  negative  residuals  were  'ob- 

(Pg-Pn)  vs  (Sg-Pg)  time  inter- 
BxLijIK:firI  with  those  for  the  earthquakes  of  2  Febru¬ 
ary  1962  and  3  March  1963  indicates  a  focal  depth  deep  in 
the  crust  for  che  former,  shallow  for  the  latter. 

OF  CONTINUOUS  VELOCITY-DEPTH  FUNCTIONS  TO 
INTERPRETATION  OF  CRUSTAL  REJECTION  DATA  IN  MISSOURI. 

S.  W.  Stewart  and  William  Stauder,  S.J.  Paper  presented  at 
Harch°19i4?al  Sooiety  of  flmerlca  meeting,  Seattle,  Wash., 

reversed  selsmic-refraotion  profiles,  each 
about  300  km  long,  were  recorded  in  northern  and  southern 
Mi  souri  in  June  1963*  as  the  result  of  a  co-operative 

fin?SueCWeo!!  ?he  U*  S*  6eol°Sical  Survey  and  St,  Louis 
lty*  Sh0t  polnts  in  northern  Missouri  reported 
rein  were  near  St.  Joseph,  Swan  Lake,  and  Hannibal.  The 
average  distance  between  recording  locations  was  11  km. 

Assuming  various  crustal  models  characterized  by  con- 

^  Jei°oltyi?Pth  Unctions,  ray  theory  was  used  to 
calculate  travel-time  curves  and  amplitudes  for  the  com- 
pressional  phases.  Several  features  of  the  observed  travel- 

^ ampllt^e  data  can  be  exPlained  by  these  ray-theory 
^alcula .ions .  For  example,  measurements  of  seismic  waves 
from  t.ie  tju?ee  shot  points  in  northern  Fils  souri  gave  three 

SH?fo0rnplltKdC'“ai?tance  curves  tor  the  Pg  phase  where  this 
phase  observed  as  a  first  arrival.  A  characteristic 1 
t*1  three  curves  is  a  sharp  decrease  of  amplitude,  fol¬ 
lowed  by  a  sharp  increase,  in  the  distance  range  140-160  km 

a  f68?}?1’6  S8*1  ?e  exPlaiRe<S  ty  changes  in  c empress ionai- 
wave  velocity  of  only  a  few  hundredths  of  a  kilometer  per 
second  in  the  upper  part  of  the  crust.  LeV~am5? litSde 

ar5i,V?JS  tbat  conHRORly  are  correctable  over  only 
8  0fi  kl^ometers  at  most  can  be  explained  by  sira- 

ilai  small  velocity  changes  within  the  earth’s  crust. 


A  REVIEW  OF  RUSSIAN  WORK  IN  MAGNITUDE  DETERMINATION 

S“nfX^o^Ch?fi^  gg;  2  *13^  D6ter- 

SeenToTe  beer, 

M  *  leg  £  -  log  £ 
rather  than 


«M  -  log  A  -  log  A* 

A  more  significant  development  is  the  wnr-ir  *hQ*.  k  ^ 
done  in  devising  an  emwv  «*oii  V  °  ,  that  h-s  been 
been  worked  out  for  the  enioent-r-'i  ?  Pra°bical  method  has 
whereby  an  energy  class  f  !  S  range  0<£<  700  km 

■>»««»•  gsn:i„s/g sl""’ 

K  -  log  E  (Joules) 

ItnclT™  t0  6Xtend  the  Saale  to  eploentral  die- 
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Abstract  (Cont'd) 


4*  Focal  niMchanlf'Li  detc  rminations .  Radiation 
patterns  of  c  and  S  waves  are  found  to  be  characteris¬ 
tic  of  given  earthquake  sequences  and,  to  some  degree 
of  given  earthquake  regions.  9 


5.  Focal  <lepth.  The  P  -p  Va  Sg-P-  method  has 
been  found  of  qualitative  help  In  studying  focal  depth 
o*  sn°cks  in  SE  Missouri.  The  transfer  function  for 
a  source  in  a  layered  medium  has  further  shown  that  at 
shallow  depths  the  source-crust  system  acts  as  a  rejec 
tion  filter  for  the  low  frequency  Dort-ion  of  the  spec¬ 
tra  of  body  waves . 


6.  Inversion  of  tody  wave  data.  The  transfer 
function  of  TW/Tu  (ratio  of  vertical  to  horizontal 
ground  motion)  has  been  calculated  as  a  function  of 
a  dimensionless  fre quern  y  for  representative  one-  and 
two-layer  crustal  models.  Spectra  of  observed  P  waves 
ai c1  then  compared  to  the  theoretical  curves  to  deter¬ 
mine  cru  *‘al  structure. 


7.  Effect  of  low  velocit  channels  on  travel 
times.  Parameters  of  a  low  velocity  channel  have 
been  systematically  varied  to  examine  the  circum¬ 
stance  under  which  the  channel  produces  or  does  not 
produce  a  shadow  zone.  ^BY  and  GNOME  data  were 
then  examined  to  study  lateral  variations  in  upper 
manule  velocities  throughout  the  central-,  eastern 
and  western  U.S,  ‘  * 


ity. 


8. 


Local  crustal  structure. 


seismicity,  pn  veloc- 
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